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The ratio of KsKs (KlKl) and KsKl production rates is calculated by considering K° — K° 
oscillation in J/tp — > K°K° decay. The theoretical uncertainty due to strong interaction in J/tp 
decay is completely canceled in the ratio, therefore, the absolute branching fractions of the CP 
violating processes of J/tp — > KsKs and KlKl can be cleanly and model-independently determined 
in case that J/tp — + KsKl decay is precisely measured. In the future r-Charm factory, the expected 
CP violating process of J/tp — > KsKs should be reached. It is important to measure J/tp to KsKs 
and KsKi, decays simultaneously, so that many systematic errors will be canceled. More precise 
measurements are suggested to examine the predicted isospin relation in J/tp — > KK decays. All 
results can be extended to decays of other vector quarkonia, <p, tp(2S) and T(1S) and so on. 
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In the Standard Model (SM) , CP violation arises from 
an irreducible weak phase in the Cabibbo-Kobayashi- 
Maskawa (CKM) quark- mixing matrix 0. CP violation 
has been established in both K and B systems. Cur- 
rently, all experimental measurements are consistent with 
the CKM picture of CP violation, and the CKM is, very 
likely, the dominant source of CP violation in low en- 
ergy flavor-changing processes 0]. However, the surpris- 
ing point is that the CKM mechanism for CP violation 
fails to account for the baryogenesis It is crucial to 
probe CP violation in various reactions, to see the corre- 
lations between different processes and probe the source 
of CP violation. 

In this Letter, we consider the possible CP asymmet- 
ric observation in J/tp — > K°K° decay, in which KsKs 
and KifCr, pairs can be formed in addition to KsK^. 
Within the SM, the possible CP violating decay processes 
of J/tp — > KsKs and KjKl are due to K° — K° oscilla- 
tion, here we assume that possible strong multiquark ef- 
fects that involve seaquarksplay no role in J/tp — > KsKs, 
Kt/Kl, and KsKl decays jj. The J/tp decays will pro- 
vide another opportunity to understand the source of CP 
violation. The amplitude for J/tp decaying to K°K° is 
(K°K°\H\ J/tp), and the K°K° pair system is in a state 
with charge parity C = — 1, which can be defined as 

\K°K°f=-i = i= [\K°)\K°) \K°)\K°)] . (1) 

Although there is weak current contribution in J/tp — > 
K°K° decay, which may not conserve charge parity, the 
K a K° pair can not be in a state with C = +1. The 
reason is that the relative orbital angular momentum of 
K°K° pair must be / = 1 because of angular momentum 
conservation. A boson-pair with I = 1 must be in an anti- 
symmetric state, the anti-symmetric state of particle- 
anti-particle pair must be in a state with C = — 1. This 
conclusion can also be illustrated by a direct calculation. 

For explicitness, let us denote the particle state \K°) 
with momentum p\ as \K°(pi)), and \K ) with momen- 



tum p2 as | K° (p2 ) ) . The general structure of the effective 
Hamitonian for J/tp — > K°K° decay is 

H = C w 9i7 M ~ &75)%c7 M (a' - 6'7 5 )c, (2) 

i—d,s 

where qi (i = d, s) and c denote the quarks d, s and 
c, Cw is the Wilson coefficient, and a, b, a', b' are the 
relevant coefficients for the vector and axial-vector cur- 
rents. Then the amplitude for J/tp decaying into a state 
\K ( Pl )K°(p 2 )} is 

(K°( Pl )K a (p 2 )\H\J/tP) 
= C w a'(K ( Pl )K (p 2 )\ J2 ©V(«-&71s)ft|0> 

i—d,s 

■(0\c^c\J/yj), (3) 

where (0|c7 M 75c| J/tp) = has been used. From the 
Lorentz structure of the above matrix element, we have 
the following decomposition 

(K ( Pl )K°(p 2 )\ £ fcy(a-frY6)?i|0) 

i— d,s 

= F + ( Pl +p 2 y + F-{ Pl - P2 r, (4) 

where F + and F~ are Lorentz invariant form factors. For 
the vector current induced vacuum- J/tp matrix element, 
there is the common decomposition 

(0|c7^c|J/V>) = fj/^mj/^, (5) 

where /ju is the decay constant, and e M the polarization 
vector of J/tp. 

Substitute Eqs.QJ and J^Jl into Eq.JSJ), and use e- (pi + 
p 2 ) = e ■ p = 0, where p is the momentum of J/tp, we can 
obtain 

{K\ Pl )K Q {p 2 )\H\J/tl:) = C w a'fj^mj^F~e-( Pl -p 2 ). 

(6) 
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Exchange p\ and p 2 , we can get another amplitude 
(K ( Pl )K°( P2 )\H\J/ij} = C w a'fj^mj^F-e-(p 2 - Pl 

The above two equations directly give 

(K ( Pl )K (p 2 ) + K a ( Pl )K (p 2 )\H\J/^) = 0, (8) 



(7) 



which shows that the amplitude for J/ip decaying into 
K°K° pair with charge parity C — +1 is zero. Therefore, 
only state with C = — 1 can be produced in J/ip — > K°K° 
decay. 

Note that Eq. 10 is based on the factorization approx- 
imation. However, the conclusion that the K°K° pair 
produced from a vector meson decay must be in C = — 1 
state does not depend on the factorization approxima- 
tion, it can be drawn directly from the momentum con- 
servation. The above computation should be only viewed 
as a complementary illustration rather than demonstra- 
tion. 

Next we shall analyze the time-evolution of K Q K° sys- 
tem produced in J/ip decay. 

The weak eigenstates of K° — K° system are \Ks) = 
p\K°) + q\K a ) and \K L ) = P \K°) - q\K°) with eigenval- 

ues fx s = ms - -Ts and PL = m L - -T L , respectively, 

where the mg and Tg (rriL and T^) are the mass and 
width of K s (K L ) meson. Following the J/ip -> K°K° 
decay, the K° and K° will go separately and the time- 
evolution of the particle states \K°(t)) and \K°(t)) are 
1 



\Ks) 



\K L )) and 



given by \K°(t)) = — [t 
2 P 

\K°(t)) = ^ {e-^ st \K s ) - e- l ^\K L )), respectively. 
Then the time-evolution of K° — K° system with C = — 1 



is 



C=-l 



lrt u (tl,i 2 )) 

[l^r°(* x )>j^°(* 2 )> - i^°(*x)>j^ (t 2 )>] 

= ^-[9ls\K l )\Ks) - gsL\K s )\K L )}, (9) 

where g LS = e^* 1 "^ 8 ' 2 and g SL = e -vsti-^t 2 . 
Since the states \Ks) and \K£) are un-orthogonal, we 
have (K S \K L ) = (K L \K S ) = \ P \ 2 - \q\ 2 and (K S \K S ) = 
(Kl\Kl) = 1. Then the amplitudes to find KsKs, 
K$Kl, KlKs and KiKl pairs are 

A 1 (t 1 ,h) = (KsKs | K°K° (ti,t2)) C= ^ 1 

' r '\p\ 2 -\q\ 2 )(9LS-gsL)}, (10) 



2V2pq 



A 2 (t u t 2 ) = {KsK L \K a K {t u t 2 )) 
1 



c=-i 



2^ P q 



[ghs 



q\ 2 ) 2 gs L ] 



(ii) 



A 3 (h,t 2 ) = (KLKslK^'it^h)) 



\C=-1 



2\f2pq 



\i\ 2 ) 2 gLs - gsh] 



1 r„ 



C=-l 



2y/2 P q 



\P\' 



q\ 2 )(gLs - 9sl)] 



(12) 



(13) 



Therefore, one can find Ai(ti,t2) = A±(t\,t2) and 
Aa(ti,t2) — —A 2 (t 2 ,ti). If CP is conserved, i.e., \q/p\ = 
1, or the two particles are observed at the same time, 
namely, t\ — t 2 , then Ai(t\,t 2 ) = An{ti,t 2 ) = for 
KsKs and KiKl cases because Bose-Einstein statistics 
prevents two identical bosons from being in an antisym- 
metric state. In other words, only Ks and Kl can be 
seen at the same time in J/ip — > K°K decay. 

Squaring the amplitudes Ai(ti,t 2 ys, one can get the 
time-dependent possibilities to find KsKs, KsKl and 
KiKl pairs 



d 2 P[if s (ti),^s(t 2 )] 



dt\dt 2 



.M s \A x {t x ^)f 



■■AT-, 



A\pq\ 2 
- cos(xT(t 2 —ti)], 

d 2 P[if L (<i),^5(i 2 )] 



1\ ) e- r ^ +t ^[cosh(yr(t 2 -h) 



(14) 



dt\dt 2 



Air 



8\pq[ 



. e -r(ti+t 2 )r e yr(t 2 -ti) 



-2(W 2 -M 2 ) 2 



x cos (xr(i 2 -t 1 )) + (H 2 -| (7 | 2 ) 4 e -^-* 1 )], (15) 



d 2 P[if S (ti),^L(^)] _ 



dt\dt% 
= N f \A 2 (t 2 M)\ 2 , 



= U f \A^t x ,t 2 )Y 



(16) 



d 2 P[X L ftx),g L ft 2 )] _ 

dt\dt 2 
A/)|A(*i,* 2 )| 2 , 



ATf\M(ti,t 2 )? 



(17) 



where TV/ is a common normalization factor, T = 

r g + r L Am Ar . 

, x = and ?/ = Am is the mass 

2 r tf 2r v 

difference of .K^, and /Cg, i.e., Am = mj, — ms, while 
Ar = — Ts is the width difference). 

The time-integrated possibilities to observe KsKs, 
KsKl and KiKl pairs, which are normalized by the 
widths of Ks and Kl, i.e., TV/ = r^r^, are 

/>oo poo 

V(K s ,Ks)=T s T L J dt 2 J *i|Ai(ti,* a )| a 

rsr^ (b| 2 - M 2 ) 2 f°i 1 y 



r 2 4|p g | 2 Vi -y 2 ! + « 
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V(K S ,K L ) 

r s r L / dt 

Jo 

T S T L 1 



3 LL I OEa / <ftl(|A 2 (tl,t2)r + p4 3 (t 1 ,t 2 )|' 

o Jo 



r 2 4| M | 2 

1 2(b| 2 



kl 2 ) 



2^2 



i-y 2 



T(K L ,K L ) 



1 + x 2 



T S T L 



f ~|g| 2 ) 4 
1-y 2 



(19) 



f/fo 



d*i|^4(ti,t 2 )| s 



= V{Ks,K s ). 



(20) 



The final amplitude squared for the decay process 
J/?/ 1 — > K1K2, where K\Ki can be KsKs, KsKl or 
KjJCl, which is from the time-evolution of K° — K° sys- 
tem, is 

\A(J/i> -> ^!^ 2 )| 2 = ViK^K^K^W/^, 

(21) 

where V(Ki,K 2 ) can be obtained from Eqs. l(T%|) . lfT^|) 
and (PU|) . Finally, the partial widths of the J/ip — > KsKs, 
KsKl and KlKl are 

rty -> X 5 ;f s ) = V(K S ,K S ) ■ C ■ \(K K°\H\^)\ 2 , (22) 



r(V> - ^s^l) = T(K S ,K L ) ■ C ■ \{K°K°\H\i>)\ 2 , (23) 



r(V» -> KlKl) = V(K L ,K L ) ■ C ■ \(K°K°\H\tP)\ 2 , (24) 



where C = — 



3 8tt m 2 /j/; 



is the phase space factor, mj/^ 



is the J/tp mass and P is the three- momentum of the 
final particles in the rest frame of J/tp- Combining 
Eqs. P0 )l . and (23 jl . one can get r(J/-0 -> Ks*f s ) = 
r(J/V -> ^l^l). The ratio R ss (R LL ) of ^ s 
(KlK l ) and K$K L production rates is obtained 



Rss 



TjJHj -> gggg) 



2\2 • 



(Ipl-lr) 



2/ 



and 



Rll 



r( J/V> -> g^i) 

r(j/^ -> ^ s ^l) 



i? 



ss- 



(25) 



(26) 



In the ratios, the phase space factor C and the strong 
matrix element squared \ (K°K°\H\J/ip}\ 2 are completely 
canceled, which ensures that the ratios are completely 
free from uncertainty caused by strong interaction in 
J/tp — ► K°K° decay. The expected value of the ratio, 
Rss (Rll), can be determined by using experimental 
measured values of x, y and \p\ 2 — \q\ 2 , where \p\ 2 — \q\ 2 



is related to the CP asymmetry parameter Sl in semilep- 
tonic Kl decay 

T(K L -> l+i 



Sl 



- T(K L -► 1-VTT+) 



T(K L 
\P\ 2 - 



l + h>n 



T(K L 



1-V1T+) 



(27) 



One can further express Rss (Rll) completely in terms 
of mixing and CP asymmetry parameters in kaon decays, 
by combining Eqs. 1 (25) 1. and l(77)l 



Rss — Rll = (Sl) 2 



1 



(28) 



where the experimental result of Sl is (3.27 ± 0.12) x 
10 -3 [jj. From the measured values of mass difference 
and widths of Kl and K$ H, one can get x — 0.946 and 
y = —0.997. The uncertainties on x and y are negligible. 

The 



Therefore. 



x 2 + y 2 



1 and R ss = Rll ~ (Sl) 



dominated uncertainty of Rss (Rll) is from the error of 
the measured value of Sl- The values of the ratios Rss 
and Rll can be obtained 



Rss = Rll = (10.66 ± 0.78) x 10" 



(29) 



which are model-independent and theoretically clean. 
The total uncertainty of Rss (Rll) is about 7%, and 
can be improved in the future if more precise measure- 
ment of CP violation parameter in Kl V^vit^ decays 
is obtained. In general, one can measure J/ip to KsKs 
and KsKl decays simultaneously in the future r-Charm 
factory, so that some of the systematic errors can be can- 
celed in the ratio. It is very interesting that the ratio can 
be extended to other quarkonia which can decay into KK 
final states, for example, in <p, tp(2S) or T(1S) — > K°K° 
decays, we have 



Rss 



r( J/y. -> KsK s ) = TW2S) -> K S K S ) 
T(J/tp -> KsK L ) T(iP(2S) -> K S K L ) 
T(</> -> KsK s ) = T(Y(1S) -> K S K S ) 

r(</> -> k s ^l) r(r(i5) -> k s ^l) 



(30) 



With current available experimental data, such as, 
B(J/V> -> ^s^i) = (1.82 ± 0.04 ± 0.13) x 10" 4 f 
B(^(2S*) -> K 5 ff L ) = (5.24 ± 0.47 ± 0.48) x 10" 5 
and B((p -> tfgK'i) = (33.7 ± 0.5)% 5], the branching 
fractions of CP violating decay processes of J/tp, tp(2S) 
and <p — * K sK s(K lK l) can be extracted as following 



B(J/V> -> i^s^s) = (1-94 ± 0.20) x 10~ 9 , 
B(ip(2S) -> K S K S ) = (0.56 ± 0.08) x 10~ 9 , 
B(4> -> ifs^s) = (3.59 ± 0.27) x 10~ 6 , 



(31) 



which arc the first model-independent predictions in rare 
quarkonium decays. With more precise measurements of 
KsKl decays, the above errors can be reduced further. 
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It is interesting that the (f> — ► KsKs decay can be even 
reached at the current KLOE experiment at the DA0NE 
accelerator. 

The CP violating decay processes of Jji\) — > KsKs 
and ip(2S) — > K$K S had been searched for by BESII 
Collaboration. The upper limits on the branching frac- 
tions at 95% C.L. are set: B(J/ip -> < 1.0 X 10" 6 
and B(ip(2S) -> #<?*Cs) < 4.6 x 1CT 6 , respectively @. 
The current bounds of the production rates are beyond 
the sensitivity for testing Rss- However, the BESIII ex- 
periment will start to take data in the middle of 2007. 
About 10 x 10 9 J/ip and 3 x 10 9 ip{2S) data samples 
can be collected per year's running according to the de- 
signed luminosity of BEPCII in Beijing 0. Thus, both 
J /if) — > KgKg and ij)(2S) — > K S K S will be reached with 
data taking in a few years at BEPCII. The cf) — > KsKs 
will also be easily accessible at the future DA0NEII Fras- 
cati (/(-factory at which the designed luminosity is about 
1.0 x 10 34 cm _2 s _1 10], and the signal is clean and free 
from backgrounds because it is just near the KK thresh- 
old. Y(1S) — > KsKs and KsKl can be studied at B- 
factory jlOj, for example, with the current luminosity at 
KEK-B, about 4.0 x 10 9 T(15) events can be collected 
with one year's running. It will be very interesting to 
collect more data at the future Super-B factory to test 
the R S s- _ 

Due to isospin symmetry, we have (K°K°\H\J/ip) — 
(K + K~ \H\J/i/j), then, it is straightforward to obtain the 
following relation by neglecting the phase space difference 



B(J/tl> -> K S K L ) 



= A ■ B(J/tp — » K + K~ 
= B{J/ip^ K+K~), 



(32) 



where A is the correction factor due to K° — K mixing, 
which can be derived from Eqs. (|19|l and (|23(1 



A 



TsT 



Si L 



1 



1 



T 2 4| M | 2L l-y 2 

2(bl 2 -ki 2 ) 2 (N 2 

1 + x 2 
0.972, 



\q\ 2 ) 4 



i 



r.sTz 



i 



(33) 



1. The 



where we have used the relation 

TV l _ yl 

uncertainty on A is tiny, at the order of 10 . The cor- 
rection factor A is slightly smaller than 1 , which is caused 
by K° — K° mixing. 

Eqs.lO and show that the ratio A = 

B(J/^ -> K S K L 



, , — ^ 1 even if isospin symmetry is an 

B{J/ip — > K+K ) 

exact symmetry. However isospin symmetry is only an 
approximate symmetry, it is usually violated by a few 
percent level. Electromagnetic effects may be even par- 
ticularly important for the the K + K~ final state. There- 
fore, the true value of the factor A may be even largely 
different from unity, which should be the total effects of 



K° — K° mixing and isospin violation. The isospin violat- 
ing effect must be taken into account when comparing the 

experimental value of the ratio A = j-f- f ^\ 

B(J/ip — * K + K~) 

with the theoretical prediction. However further study 
of the isospin violating effect is beyond the scope of this 
paper. 

One more remark is in the following. Soft photons 
can be emitted from the initial and final states in vec- 
tor quarkonia to KK decays. The radiation of the soft 
photons in the decays allows the K°K° in a C = +1 
state. Such a process with a soft photon in the KsKs 
or K[Kl final state is not CP-violating. The detection 
of the soft photons depends on the sensitivity of the de- 
tectors. Therefore the soft-photon-radiation process is 
an experimental background for the test of the Rss an d 
Rll predictions. This background should be subtracted 
in experiment. 

In conclusion, we have studied the CP violating decay 
processes of J/tp, tp(2S), <j> and T(1S) quarkonia to KsKs 
and KlKl- The ratio R SS (Rll) of K S K S (KlK l ) 
and KsKl production rates has been constructed in a 
model-independent and theoretically clean way. Simul- 
taneous measurements of vector quarkonium decays to 
both KsKs and KsKl pairs are suggested at higher lu- 
minosity e + e~ machines, so that many systematic errors 
can be canceled. With the current experimental informa- 
tion, the absolute branching fractions of the CP violating 
processes are firstly predicted. The isospin relation is ob- 
tained by considering K° — K° mixing effect. 
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